Structurally-driven magnetic state transition of biatomic Fe chains on Ir(OOl) 
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Using first-principles calculations, we demonstrate that the magnetic exchange interaction and 
the magnetocrystalline anisotropy of biatomic Fe chains grown in the trenches of the (5 x 1) re- 
constructed Ir(001) surface depend sensitively on the atomic arrangement of the Fe atoms. Two 
structural configurations have been considered which are suggested from recent experiments. They 
differ by the local symmetry and the spacing between the two strands of the biatomic Fe chain. 
Since both configurations are very close in total energy they may coexist in experiment. We have 
investigated collinear ferro- and antiferromagnetic solutions as well as a collinear state with two mo- 
ments in one direction and one in the opposite direction (|J.|-state). For the structure with a small 
interchain spacing, there is a strong exchange interaction between the strands and the ferromagnetic 
state is energetically favorable. In the structure with larger spacing, the two strands are magneti- 
cally nearly decoupled and exhibit antiferromagnetic order along the chain. In both cases, due to 
hybridization with the Ir substrate the exchange interaction along the chain axis is relatively small 
compared to freestanding biatomic iron chains. The easy magnetization axis of the Fe chains also 
switches with the structural configuration and is out-of-plane for the ferromagnetic chains with small 
spacing and along the chain axis for the antiferromagnetic chains with large spacing between the 
two strands. Calculated scanning tunneling microscopy images and spectra suggest the possibility 
to experimentally distinguish between the two structural and magnetic configurations. 



I. INTRODUCTION 

Driven by the wish to realiz e the proposed concepts of 
future spintronic devicea^^ the development of novel 
nanostructures and nanomaterials with tailored elec- 
tronic and magnetic properties has become a key chal- 
lenge of today's research. A promising path to control the 
magnetic properties of matter is to use low-dimensional 
systems and to reduce their size down to the nanome- 
ter or even atomic scale. For nanoscale systems, how- 
ever, an essential requirement is to enhance the magnetic 
anisotropy in order to stabilize magnetic order against 
thermal fluctuations or quantum tunneling. The manip- 
ulation of exchange interactions opens another path to 
create new materials with a magnetic state that may be 
tunable by external magnetic or electric fields!^ E.g. the 
occurrence of chir al sp in spiral states at surfaces has 
been demonstrated, 5 6 their manipulation by electrical 
currents was suggested 5 -*^ and a way to grow films with 
multiple metastable magnetic states has been proposed P 

The ability to create one-dimensional monoatomic 
magnetic chains of transition-metals on surfaces by self- 
organizatiorpiifil or by manipulation with a scanning tun- 
neling microscope^ has recently opened new vistas to 
explore and manipulate artificial magnetic nanostruc- 
tures even atom-b y-ato m. E.g. the pioneering work of 
Gambardella et a/P^ demonstrated that the magnetic 
anisotropy of atomic transition-metal chains, consist- 
ing of Co atoms on a stepped Pt(lll) surface, is dra- 
matically enhanced with decreasing dimensionality from 
two-dimensional films to quasi-one-dimensional chains 
and depends sensitively on the number of Co strands 



in a chain. These experimental observations were ex- 
plained based on electronic structure calculations which 
emphasized the crucial role played by the substrate, 
reduced symmetry, and struc tural rela xations for the 
magneto-crystalline anisotropy! 12 * 13 * 1 4 ^ The large mag- 
netic anisotropics led to slow relaxation dynamics of the 
magnetization and the observation of magnetic hystere- 
sis loops at low temperatures indicative of ferromagnetic 
coupling. In another experiment, Mn chains of up to 10 
atoms were created by manipulation with a scanning tun- 
neling microscopy tip on an insulating CuN layer grown 
on Cu(001). Experimentally, the exchange interaction 
between individual spins was obtained by measuring the 
excitation spectrum via inelastic tunneling spectroscopy 
which showed the quantum behavior of the entire chainP^ 
Even the sign and size of the exchange interaction be- 
tween the Mn atoms could be extracted from the exper- 
imental data. Calculations based on density-functional 
theory (DFT) clarified that a superexchange mechanism 
along the Mn-N-Mn bon d is re sponsible for the weak an- 
tiferromagnetic coupling! 16 * 17 ^ 

Recently, Hammer et al. have used a combination 
of IV-LEED and scanning tunneling microscopy (STM) 
measurements to demonstrate that the Ir(001) surface 
can serve as an ideal template to grow defect-free, 
nanometer long transition-metal nanowires of different 
structure, chemical composition, and le ngth depend- 
ing on the preparation conditions! 18 * 19 * 21 ^ E.g. biatomic 
Fe chains can be created on the (5 x 1) reconstructed 
Ir(001) surface and lifting the surface reconstruction by 
hydrogen opens the possibility to produce Fe-Ir-Fe tri- 
atomic chains. While the biatomic Fe chains are a very 
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promising system to study magnetism of (quasi-) one- 
dimensional transition-metal chains, there is little un- 
derstanding so far. Experimentally, it is extremely chal- 
lenging as for laterally averaging measurements samples 
with a homogeneous distribution of chains are needed 
or a technique must be applied which allows to locally 
probe the magnetic properties of individual Fe chains. 
Another key difficulty is the detailed characterization of 
the chains structure. From combined STM and LEED ex- 
periments it is only known that the biatomic chains grow 
in the trenches of the (5 x 1) reconstructed Ir(OOl) sur- 
face, however, the adsorption sites in the trenches could 
not be deduced. As structure and magnetism are closely 
correlated in such systems, their magnetic properties are 
an open issue. 

A theoretical study using first-principles calculations 21 
reported an excellent agreement with the structural pa- 
rameters of the (5 x 1) reconstructed Ir(OOl) surface and 
concluded that the Fe chains are strongly ferromagnetic 
at low temperatures but were probably non-magnetic 
in the room temperature measurements of Hammer et 
al. However, this theoretical study considered only fer- 
romagnetic solutions and did not determine the magne- 
tocrystalline anisotropy energy which is crucial for an ex- 
perimental verification of the proposed ferromagnetism in 
these chains. On the other hand, a bd transition-metal 
substrate such as Ir can have dramatic consequences on 
the exchange coupling in a deposited Fe nanostructure as 
is apparent from the observation of a complex nanoscale 
magnetic structure for an Fe monolayer on Ir(lll)j22J 
the antiferromagnetic ground state of an Fe monolayer 
on W(001), 23 to name just a few. Recent studies on 
Fe stripes on Pt(997)P^ and FePt surface alloys^ report 
on strong correlation between complex magnetic ground 
states and the details of structural arrangement. 

Here, we use first-principles calculations based on 
density-functional theory to study the structural, elec- 
tronic, and magnetic properties of biatomic Fe chains 
deposited in the trenches of the (5 x 1) reconstructed 
Ir(OOl) surfacep3 We focus on two structural arrange- 
ments of the biatomic chains which differ by the ad- 
sorption sites of the Fe atoms. In one configuration the 
distance between the two strands of the biatomic chain 
is smaller than the atom spacing along the chain direc- 
tion (denoted as CI in accordance with Ref. |2"T|| . while 
in the other their separation is clearly larger than the 
interchain spacing (denoted as C4), see Fig.[l] We con- 
sider collinear ferro- and antiferromagnetic arrangements 
along the chains and the | It-state with two moments in 
one and one moment in the opposite direction. We find 
that the energetically favorable magnetic state as well as 
the easy magnetization axis of the Fe chains depend sen- 
sitively on their atomic arrangement and local symmetry. 
For the Fe chains in the Cl-structure, the exchange cou- 
pling along the chain is ferromagnetic, however, due to 
hybridization with the Ir substrate it is much weaker than 
for freestanding biatomic Fe chains with the same atom 
spacing. Due to their small separation, the two strands 



of the Cl-chains are also ferromagnetically coupled. Sur- 
prisingly, for the C4-configuration, we find a transition 
from ferromagnetic state for free-standing chains to anti- 
ferromagnetic order along the biatomic chains after depo- 
sition on the Ir substrate. In this case, the hybridization 
with the Ir surface is strong enough to invert the sign of 
exchange coupling, while the two Fe strands are nearly 
exchange decoupled. The interplay of the Fe interstrand 
distance and the hybridization with the Ir substrate re- 
sults also in a different easy axis of the magnetization 
for the two structures: while in the Cl-FM state the 
easy axis is out-of-plane, it switches into the chain axis 
for the C4-AFM configuration. The total energy of the 
CI- and C4-structure are quite close and therefore both 
chain types could occur in an experiment depending on 
the growth conditions. We simulate measurements by 
STM and observe that the two strands of the biatomic 
Fe chain in the Cl-ferromagnetic state are too close to be 
individually resolved, while they can be distinguished in 
the C4-structure with an antiferromagnetic ground state. 
In the latter case, spin-polarized STM (SP-STM) should 
further allow to directly resolve the two-fold magnetic 
periodicity along the chain. 

The structure of this paper is as follows. In the next 
section we present details of our method and the calcu- 
lations. Then we discuss the structural relaxations of 
the pure (5 x 1) reconstructed Ir(OOl) substrate and of 
the biatomic Fe chains in the two structural configura- 
tions on the (5 x 1) reconstructed Ir(OOl) surface. In sec- 
tion IV, we analyze the magnetic ground state configura- 
tion and the effects of hybridization with the Ir substrate, 
before we turn to the magnetocrystalline anisotropy in 
section V. We explore the feasibility to experimentally 
resolve the different structural and magnetic properties 
by SP-STM and to verify our predictions of a structure- 
dependent magnetic ground state. Finally, a conclusion 
and summary is given. 



II. COMPUTATIONAL DETAILS 

We employed the film- version of the full-potential lin- 
earized augmented plane-wave (FLAPW) method, as im- 
plemented in the Jiilich density-functional theory (DFT) 
code FLEUR. We used inversion-symmetric films with 7 
and 9 layers of the (5 x 1) reconstructed Ir(OOl) sur- 
face and Fe biatomic chains on both sides of the slab. 
The whole system possesses spatial inversion symme- 
try. We calculate the biatomic Fe chains in a (5 x 1) 
supercell along the y-direction ([011]-axis), perpendicu- 
lar to the chain x-axis ([011]-axis), which results in a 
distance of 13.75 A between the axes of two adjacent 
biatomic chains, c.f. Fig. [T] for structural arrangements 
and definition of the axes. The theoretical Ir lattice 
constant we used for calculations constituted a value of 
3.89 A. We used the generalized-gradient approximation 
(GGA, revPBE functional^) of the exchange-correlation 
potential for the structural relaxations and tested total 





FIG. 1: (color online) Geometrical structure of the biatomic Fe chains on the (5 x 1) reconstructed Ir(OOl) substrate in CI 
((a)-top view, (c)-side view) and C4 ((b)-top view, (d)-side view) configurations. Fe atoms are marked in blue, while Ir atoms 
are marked in gold. In (e) and (f) a schematic top view of respectively CI and C4 biatomic Fe chains is shown (positions of 
the atoms do not correspond to realistic calculated values). In the latter graphs the numbering of the Ir atoms corresponds to 
that in Fig. 2 and table II. The x- and j/-axis defined in (a) correspond to [Oil] and [Oil] directions, respectively. 



energy differences between the two magnetic configura- 
tions also within the local-density approximation (LDA, 
VWN functional^) . We used 18 fc-points in a quarter 
of the full two-dimensional Brillouin zone (2D-BZ) for 
self-consistent calculations. The calculated total energy 
differences between different magnetic ground states were 
carefully tested with respect to the number of fc-points. 
For the basis functions, we used a cut-off parameter of 
/: max = 3.6 a.u." 1 for relaxations and 3.7 a.u. -1 for com- 
paring the total energies of different magnetic configura- 
tions. 

We considered two possible structural arrangements of 
the Fe atoms denoted as CI and C4 according to the no- 
tation of Ref. [HJ which are shown in Fig. [T] Experimen- 
tally, it has been observed by STIvP^ that the biatomic 
chains grow in the trenches of the (5x1) reconstructed 
Ir(OOl) surface and we therefore focus on these two con- 
figurations. Relaxations were performed until the forces 



changed by less than 3 • 1CP 4 htr/a.u. The convergence 
of the relaxed atomic positions was carefully tested with 
respect to the computational parameters. Relaxations 
were performed only for the ferromagnetic state of both 
CI and C4 configurations and the antifcrromagnetic and 
the t|t-states were calculated on these atomic positions. 



III. STRUCTURE AND RELAXATIONS 

As a first step, we performed a structural relaxation of 
the pure (5 x 1) reconstructed Ir(OOl) surface. As can be 
seen in tablc|TJ the results we obtain with a film of 9 layers 
agree very well with the experimental data measured by 
IV-LEEDp3 In particular, the experimentally observed 
trench-like structure is reproduced. Our values are also 
in close agreement with those obtained with the VASP 
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FIG. 2: (color online) Schematical representation of the 
structural parameters given in table I. Note that in this figure 
the buckling of the Ir substrate is greatly exaggerated and the 
positions of the Fe atoms (blue) and Ir atoms in the first and 
second layer of the substrate (orange) do not correspond to 
realistic calculated values. 



code by Spisak et al. 21 For deposited biatomic Fe chains 
we have also found a very good agreement between the 
values of the relative atomic positions obtained by relax- 
ing a slab with 9 and 7 layers of the Ir substrate (see 
table 0. 

We performed relaxations of the Fe biatomic chains 
only for the ferromagnetic solutions to reduce the com- 
putational effort, presenting the results in table [TJ While 
the distance A between the atomic strands in Cl- 
configuration is 2.35 A, and thus smaller than along the 
chain (see Figs. 1 and 2), it is almost twice larger in the 
C4-configuration, i.e. A = 4.23 A, and the two strands 
are well separated. In the latter structure, the Fe atoms 
are also more embedded into the Ir surface which is il- 
lustrated by the smaller vertical distance 5 from the Ir 
surface atoms between the two Fe chain atoms (denoted 
as 1x4 in Fig. [2J. The influence of the C4-chains on the 
buckling of the Ir substrate is also more pronounced as 
can be seen from the increased vertical separations of the 
surface Ir atoms. The different structural relaxations of 
the two chain configurations already hint at a larger hy- 
bridization of the Fe 3d and Ir 5d states and a stronger 
influence in the C4 arrangement . For both configurations 
the distance between the Fe dimers along the chain's axis 
was imposed by the Ir substrate and constituted 2.75 A. 



IV. MAGNETIC ORDER AND EXCHANGE 
INTERACTIONS 

Now we study the magnetic order and exchange inter- 
action of the chains in the two structural arrangements. 
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TABLE I: Relaxations of the Fe biatomic chains and the 
uppermost layer of the (5 x 1) reconstructed Ir(001) substrate 
for 9- and 7-layer slabs. All values are given in Angstrom. The 
distances d, b, and p correspond to those in Ref. 29 and are 
depicted in Fig. [2] For comparison the relaxations of the bare 
Ir(001) substrate from Ref. [29] (IV-LEED data) and Ref. ED 
are given. For the interchain distance A values is brackets 
correspond to those calculated in Ref. 1211 



We considered the ferro- (FM) and antiferromagnetic 
(AFM) solution (with antiparallel magnetic moments be- 
tween adjacent Fe atoms along the chain) for both types 
of chains. For comparison, we have also calculated free- 
standing biatomic Fe chains with the same interatomic 
spacings as in the C4 and CI configurations. As can be 
seen in table [TTJ the magnetic moments of the Fe atoms 
arc slightly larger in the C4 configuration where the Fe 
atoms are further apart. The difference in the magnetic 
moments between FM and AFM state is very small. The 
induced moments in the Ir surface layer depend much 
more sensitively on the magnetic state of the Fe chains. 
For the FM solutions, these moments are significant in 
both structures and decay slightly faster for the CI- than 
for the C4-configuration. In the AFM state, due to sym- 
metry some of the Ir atoms do not carry an induced mo- 
ment. In the C4 structure, the Ir surface atom between 
the Fe atoms has a rather large moment of 0.18 /is indi- 
cating a strong hybridization. 

The total energy differences between the two mag- 
netic configurations given in table |TT] reveal a surprising 
result. In particular, we find that the preferred mag- 
netic state is FM for the Fe chains in the CI struc- 
ture (dp c _p c —2.35 A), while the total energy difference 
is in favor of the AFM solution for the C4 structure 
with Fe atoms further apart (<iFc-Fc=4.23 A). For free- 
standing biatomic Fe wires with the same spacing be- 
tween Fe atoms, the energy difference between FM and 
AFM states is 74 meV/Fe and 164 meV/Fe in favor of 
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the FM state for the CI and C4 configuration, respec- 
tively. For supported chains the corresponding energy 
differences are of the order of rather small 30 meV/Fe, 
indicating a strong influence of the Ir substrate on the 
magnetic coupling and weakening of the FM interaction 
in the free-standing chains due to hybridization with the 
Ir atoms.^1 Because of the large separation of the two 
strands of the biatomic chain in the C4 configuration, 
this effect is dramatic and leads to the AFM ground 
state. This notion is further supported by the magni- 
tude of the exchange interaction between the two strands. 
We have calculated the total energy difference between 
a FM and AFM alignment of the two strands of the de- 
posited Fe chains which is 166 meV/Fe and 4 meV/Fe 
in favor of ferromagnetic coupling in the CI- and C4- 
configuration, respectively. The extremely small value in 
the C4-configuration, indicating that the two strands are 
magnetically nearly decoupled, is probably due to a small 
indirect exchange interaction via the Ir substrate. 

We have also calculated the FM and AFM total en- 
ergy differences for both CI and C4 configurations using 
7 layers of Ir substrate, and the results are compared 
in table |TT] to those obtained with 9 layers. In the Cl- 
structure the FM-AFM energy difference of 20.6 meV/Fe 
is in good agreement with the value of 21.1 meV/Fe for 9- 
layers of substrate. In the C4-configuration, on the other 
hand, the energy difference between the AFM ground 
state and the FM state increases by only 7 meV/Fe for 
the thinner substrate. These results show that the favor- 
able magnetic state within each structural arrangement 
is independent of the substrate thickness. However, in 
case of 9 layers the C4-AFM state is lower in energy 
than the Cl-FM solution by 9.8 meV/Fe, while in case 
of 7 layers this energy difference reverses sign and consti- 
tutes 8.5 meV/Fe in favor of Cl-FM state (see tabic [n}. 
This discrepancy probably arises due to quantum well 
states in the Ir substrate which do not influence the to- 
tal energy difference between different magnetic states 
within the same structural arrangement. Overall, our 
calculations reveal that, judging only from this energy 
difference, both configurations might appear in experi- 
ment and can be observed via, e.g., scanning tunneling 
microscopy measurements. 

We have further checked the influence of the exchange- 
correlation potential on the FM-AFM energy differences 
and found very similar results within the local-density 
approximation LDA. Using the LDA, a 7 layer (5x1) 
Ir(001) substrate and the relaxed atomic positions found 
with GGA, in the CI configuration the FM state is by 
25.1 meV/Fe lower in energy than the AFM state (c.f. a 
value of 20.6 meV/Fe in GGA). In the C4 structure the 
AFM state is by 23.7 meV/Fe lower than the FM state 
(c.f. a value of 40.3 meV/Fe in GGA). 

Within GGA, we also performed calculations of a 
collincar magnetic state with a spin arrangement of tit 
(periodically repeated along the chain axis) for the Fe 
biatomic chains on 7 layers of Ir substrate in the CI and 
C4 structure. Our calculations reveal that the C4-||t 



state is by 16.0 meV/Fe higher in energy than the C4- 
AFM state and by 24.3 meV/Fe lower in energy than 
the C4-FM state. In the CI configuration, the Tit-state 
is 9.4 meV/Fe higher in energy than the FM state, and 
by 11.2 meV/Fe lower in energy than the AFM state. 
From these three collinear magnetic solutions we can es- 
timate the nearest-neighbor and next-nearest neighbor 
exchange constants of an effective ID Heisenberg model 
to be J i = —10 meV and J2 = +1 meV for the C4 struc- 
ture and Ji = +5 meV and J2 = —1.3 meV for the CI 
structure. These values illustrate the strong tendency to- 
wards antiferromagnetic coupling due to the Ir substrate. 
Even in the Cl-configuration, the FM nearest-neighbor 
coupling has become very weak. Note, that a similar in- 
fluence of the Ir surface has recently been reported for Fe 
monolayers on Ir(lll) and for other Ad/5d substrates.^ 

In order to understand the sensitive dependence of the 
magnetic coupling in the Fe chains upon the structural 
arrangement we take a look at the density of states (DOS) 
for the two configurations, shown in Fig. [3] in compari- 
son with the free-standing chains. In the non-magnetic 
state, the DOS of the supported Fe chains displays a large 
peak at the Fermi energy for both structures, however, 
the d-band width is smaller in the C4 structure due to 
the larger separation and weaker hybridization between 
the Fe atoms — an effect even more pronounced in free- 
standing wires. 

In the CI structure, the direct interaction between the 
Fe atoms perpendicular to the chain is much stronger 
than in C4-wires and the free-standing chains are a better 
approximation. Correspondingly, the FM DOS of the 
free-standing chains in Cl-geometry is very similar to the 
supported chains, while larger changes are visible in the 
C4-FM configuration. This becomes even more evident 
from the comparison of the bandstructures of the Cl-FM 
free-standing and Cl-FM supported biatomic Fe chains 
presented in Fig. [3] In this plot the electronic states of 
the free-standing CI chain (small red and small green 
circles) display a close correspondence to the states of 
the CI supported chain (large black circles) which are 
localized mainly inside the muffin-tin spheres of Fe atoms. 
Remarkably, for many of the bands of the two systems a 
direct correspondence in terms of symmetry can be made. 

For the AFM solution in both structural configura- 
tions the modifications in the DOS due to interaction 
with the substrate are quite significant. The electronic 
bands in free-standing AFM chains are normally ver y fl at 
and corresponding peaks in the DOS are very sharp^ — 
in this case the effect of the hybridization of the local- 
ized 3d-orbitals with extended states of Ir atoms on the 
DOS can be very strong. For both FM and AFM mag- 
netic states, a slightly larger exchange splitting can be 
observed in C4 deposited chains, as compared to the CI 
configuration, which leads also to larger spin moments of 
Fe atoms in C4 arrangement (c.f. table |TT]). Overall, an in- 
terplay of decreasing Fe-Fe hybridization with increasing 
Fe-Ir hybridization when going from CI to C4 structural 
configuration leads to somewhat larger localization of Fe 
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TABLE II: Relative total energies obtained in GGA (in 
meV/Fe-atom), spin moments in the muffin-tin spheres of the 
Fe atoms in unsupported (UBC) and supported bichains, as 
well as for Ir surface atoms, and total moments in the unit 
cell (in (j,b) for calculations with 9 and 7 layers of the (5x1) 
reconstructed Ir(001) substrate. The Ir surface atoms are de- 
noted as in Fig. [2] 

electronic states in C4 biatomic chains. 



V. MAGNETO-CRYSTALLINE ANISOTROPY 

As mentioned in the introduction, the magneto- 
crystalline anisotropy energy (MAE) is a key quantity 
for nanoscale magnets as it determines the preferred di- 
rection of the magnetic moments and is crucial to sta- 
bilize the magnetic order against thermal fluctuations. 
We have calculated the MAE for the Fe chains within 
the GGA employing the force theorem and found that 
its value is stable against the chosen ground state upon 
which we perform the perturbation. For the CI structure, 
we considered the FM ground state and started from a 
converged ground state with an out-of-plane magneti- 
zation from which the MAE was obtained by applying 
the force theorem for three possible high-symmetry di- 
rections: perpendicular to the surface and the chain axis 
(along the z-axis), parallel to the axis of the chain (||- 
direction) and in the surface plane perpendicular to the 
chain's axis (Indirection) . We define two principal ener- 
gies: MAEy = E tot (||) - E tot (z), and MAEj_ = E tot (l_) 
— Etot(-z)- The number of fc-points used for calculations 
constituted 144 in the full 2D-BZ. We carefully tested 
the MAE values with respect to the number of fc-points. 
Using 144 /c-points results in an accuracy of not less 
than 0.3 meV/Fe. We find that in the Cl-FM config- 
uration the calculated values are 2.1 meV/Fe for MAEp 
and 1.8 meV/Fe for MAEj_ which corresponds to a mag- 



netization along the z-axis in the ground state, i.e. per- 
pendicular to the surface and chain axis, see Fig. [5]ja). 

For the calculations of the MAE in the C4-AFM con- 
figuration we used a slab of 7 layers of the (5 x 1) Ir(001) 
substrate with 8 Fe chain atoms in the unit cell and 144 
fc-points in the full 2D-BZ (again the stability of the MAE 
with respect to the number of /c-points was carefully 
checked). We obtain values of —0.2 meV/Fe for MAEii 
and 1.2 meV/Fe for MAEj^ corresponding to a ground 
state with the magnetic moments along the chain axis, 
see Fig. [5jb). Considering the accuracy of our calcula- 
tions we conclude that the easy axis of the magnetization 
for the Fe biatomic chains in the C4-AFM configuration 
constitutes a plane which cuts through the chain axis and 
is perpendicular to the substrate. 

The MAE and its dependence on the configuration can 
be qualitatively related to the anisotropy of the orbital 
moments in the system.^ For this purpose for differ- 
ent magnetization directions we compare the orbital mo- 
ments inside the atomic spheres of Fe chain atoms, /^ e , 
of the Ir atoms in the surface layer, /z^, and the total 
orbital moment per unit cell, /i^ ot , defined as the sum 
over the moments of the two Fe chain atoms and the five 
Ir surface atoms on one side of the slab. For the evalua- 
tion of ^£ 0t we do not take into account the much smaller 
contributions from Ir atoms deeper inside the slab. For 
the antifcrromagnetic ground state in the C4 geometry, 
one should note that the sign of the orbital moments 
switches for atoms of antiparallel magnetization. There- 
fore, we define /i^' in this case as the orbital moment 
summed over atoms in one half of the (2 x 5) unit cell on 
one side of the slab. 

In the Cl-FM state the anisotropy of fj,^ is in quali- 
tative agreement with the anisotropy of the total energy: 
while for an out-of-plane magnetization /i^ 3 * reaches 
a value of 0.126/^b, it constitutes only 0.083/is and 
0.035/is for the _L- and ||-directions of the magnetiza- 
tion, respectively. Therefore, the easy axis coincides with 
that along which the orbital moment is largest p3 This 
anisotropy of ^ % can be explained based on the depen- 
dence of the Ir contributions on the magnetization di- 
rection. For an out-of-plane magnetization, the values 
of fji^l are small and the total orbital moment is domi- 
nated by the Fe atoms. In contrast, for a magnetization 
along the _L- and the Indirection the Ir orbital moments, 
H l £, reach significant values, however, of opposite sign 
with respect to /i^* 3 . For the ^-magnetization direction 
with /z£ e = 0.077/1 s the value of the orbital moment 
of the 1x4 atom (see Fig. 2) is — 0.013/is, while for the 
1 1 -magnetization it even reaches — 0.029/is and the cor- 
responding value of is only 0.058/x_b- 

In the C4-AFM configuration of the Fe chains the 
agreement between the anisotropy of the orbital moment 
and of the total energy is even better. While for the 
^-magnetization direction the value of is 0.160/is, 
it reaches much larger values of 0.239^b and 0.253/is 
for the z- and || -direction, respectively. The direction of 
the smallest total orbital moment, /i^ ot , coincides with 
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FIG. 3: (color online) Density of states (DOS) for the nonmagnetic (NM), ferromagnetic (FM) and antiferromagnetic (AFM) 
state of Fe biatomic chains on 9 layers of the (5 x 1) reconstructed Ir(001) surface in C4, (a)-(c), and CI, (d)-(f), configuration. 
The DOS is given inside the muffin-tin spheres of Fe atoms. The dashed lines show the DOS for unsupported Fe biatomic 
chains with corresponding spacings between the Fe atoms in the chain. In (b), (c), (e) and (f) left (green) and right (red) 
curves stand for spin-up and spin-down channels, respectively. 



the hard axis and, moreover, a very small difference be- 
tween the values of /i^ ot for the two other directions cor- 
responds to a very small energy difference MAE|| . In this 
structural arrangement the anisotropy of the Fe orbital 
moments, /x^ e , dominates the anisotropy of /i^ ot . The 
contribution of the Ir atoms to the total orbital moment 
is nearly independent of the magnetization direction. For 
the z- and || -magnetization directions /x£ e is 0.100/is and 
0.106/is, respectively, while it is only 0.060/is and thus 
much smaller for the hard _L-magnetization axis. 



VI. SIMULATION OF STM EXPERIMENTS 

In the previous sections, we have demonstrated that 
the magnetic properties of biatomic Fe chains on the 
(5 x 1) Ir(OOl) surface depend crucially on the atomic 
arrangement of the atoms. Both easy magnetization di- 
rection and magnetic order change upon displacements 
of the Fe atoms which changes the interaction between 
the two strands of the chain and their hybridization with 
the Ir substrate. In order to verify our predictions ex- 
perimentally, a technique with a high lateral resolution 
seems indispensable. Therefore, we study theoretically 



the possibility to use spin-polarized scanning tunneling 
microscopy (SP-STM) to resolve the atomic and mag- 
netic structure of these chains. 

Fig. [6] displays the calculated local density of states 
(LDOS) in the vacuum at a distance of about 6.8 A from 
the Fe chain atoms. Within the Tersoff-Hamann model 
of STM the vacuum LDOS is directly comparable with 
measured dl/dll tunnel spectra. The comparison of dif- 
ferent magnetic states shows distinct features in the FM 
solution for both chain structures. Strong peaks appear 
in the minority spin channel at +0.7 eV and +0.4 eV in 
the C4 and CI configuration, respectively, see Figs. |6jc) 
and (d). In the AFM state, we observe a structure of 
two broad peaks at —0.5 eV and +1.0 eV for the CI 
chains, Fig. [6^f ) , while the C4 chains display a relatively 
featureless structure with a small peak at the Fermi en- 
ergy. Considering the two favorable Cl-FM and C4-AFM 
states of the chains, it seems that the two types can be 
distinguished by the strong peak of the Cl-FM configu- 
ration. 

A direct way of verifying the magnetic structure of the 
Fe chains might be feasible by imaging them in the to- 
pography mode of SP-STM. In addition, the structural 
arrangement of the atoms might be detectable. Our sim- 





X k r k X 

FIG. 4: (color online) Bandstructure of the Fe biatomic 
chains in the Cl-FM configuration on 9 layers of the (5 x 1) re- 
constructed Ir(OOl) surface. Left panel shows majority bands 
and right panel the minority bands, k- vector has been cho- 
sen along the chain direction. Large black filled circles denote 
states which are localized on the Fe chain. For comparison the 
bands for an unsupported biatomic Fe chain with the same 
interatomic spacing are given by green (majority) and red 
(minority) filled small circles. Grey circles on the background 
mark the states of the whole Fe+Ir system. 



ulations of STM and SP-STM images are displayed in 
Figs. [7] and [8] for the FM ground state of the CI chains 
and the AFM ground state of the C4 chains, respectively. 
We have chosen an energy window corresponding to the 
unoccupied states close to the Fermi energy, but the oc- 
cupied states lead to very similar results. 

The two strands of the biatomic Fe chains in the CI 
structure are only 2.35 A apart which makes it impos- 
sible to resolve them in cross-sectional scans as seen in 
Fig. [7]ja) and (b) for both spin channels and also in their 
summation. The corrugation amplitudes which we calcu- 
late from these plots are 1.9 and 2.5 A and the apparent 
width (full width at half maximum using the topmost 
lines of constant charge density) amounts to 7.7 and 6.9 
A for the majority and minority spin channel, respec- 



FIG. 5: (color online) Sketches of the energy landscape as 
a function of magnetization direction for the two types of 
chains, (a) Cl-FM configuration with an easy axis pointing 
perpendicular to the surface, hard axis along the chain (+2.1 
meV/Fe atom), and middle axis perpendicular to the chain 
and in the surface plane (+1.8 meV/Fe atom), (b) C4-AFM 
configuration with an easy axis along the chain axis but only 
a small energy difference with respect to the middle axis per- 
pendicular to the surface (+0.2 meV/Fe atom), and a hard 
axis perpendicular to the chain and in the surface plane (+1.2 
meV/Fe atom). In this plot, we construct the energy land- 
scape based on the lowest power of the directional cosines of 
the magnetization with respect to the crystallographic axes 
which are allowed by symmetry. The coefficients are obtained 
from the values of MAEy and MAE±. 

tively. Interestingly, the chains appear wider if major- 
ity states of the chain are imaged than for the minority 
states which can be explained based on the orbital char- 
acter of the dominating states. From the DOS of the 
Fe atoms, Fig. [3je), we see that the minority channel 
is dominated by d-electrons while the majority d-band 
is far below the Fermi energy and therefore, s- and p- 
states provide a large contribution. The (i-character of 
the states for spin down electrons is clearly visible in the 
cross-sectional plot, Fig.^b), and leads to a sharper im- 
age of the chains. The more delocalized s- and p-states 
dominate the majority spin channel, Fig. [7ja), and lead 
to a larger apparent width and smaller corrugation am- 
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FIG. 6: (color online) Local density of states (LDOS) in 
the vacuum for the NM, FM and AFM states of Fe biatomic 
chains on 9 layers of the (5x1) reconstructed Ir(001) surface in 
C4 (left column) and CI (right column) configurations. The 
LDOS is given at a distance of ~ 6.8 A above the Fe atoms. 

plitude of the chains. 

In the C4 configuration, the Fe chain atoms are much 
further apart in the perpendicular direction, d = 4.23 A, 
which is large enough to allow the resolution of the two 
strands as can be seen in our simulations of STM images, 
shown in Fig. [8] The corrugation amplitude of the C4- 
AFM chains in both spin channels is 1.5 and 2.2 A in the 
spin up and down channel, respectively, the corrugation 
between the two Fe chain atoms is 0.15 and 0.3 A, and 
the apparent widths are about 10.5 A, i.e. significantly 
wider than the CI FM chains. With a spin-polarized 
STM it should further be possible to resolve the antiferro- 
magnetic spin alignment along the chains as can be seen 
from the STM images for the two separate spin channels 
as seen from Figs. [8^c) and (d). For a spin-polarization 
of the tip of Pt — (nj — njj/(n.f = 0.4, where n-\ and 
ni are the majority and minority spin LDOS of the tip 
at the Fermi energy, we obtain a corrugation amplitude 
of Az = 0.15 A along the chain, i.e. the maximum height 
change as the tip scans along the chain. 

VII. CONCLUSIONS 

We performed extensive first-principles calculations to 
elucidate the interplay of structure and magnetism in bi- 



atomic Fe chains on the (5 x 1) reconstructed Ir(001) 
surface. We find a crucial influence of the hybridiza- 
tion of Fe chains with the Ir substrate on the magnetic 
ground state of the wires. Depending on the particu- 
lar structural arrangement, the magnetic ground state 
switches from along-the-chain ferromagnetic for the CI 
configuration with a smaller (w 2.4 A) distance between 
the two strands, to antiferromagnetic for the C4 state 
for which this distance constitutes an almost twice larger 
value (w 4.2 A). In the C4 configuration, the two strands 
of the chain are nearly decoupled in terms of exchange 
interaction, while we find strong ferromagnetic coupling 
in the CI configuration. We also find that the direc- 
tion of the magnetization in these two configurations is 
different: while in Cl-FM chains the Fe spin moments 
point out of plane with a value of the magnetic anisotropy 
of ~ 2 meV/Fe with respect to in-plane directions, the 
magnetization in C4-AFM chains can freely rotate in the 
plane of along-the-chain and out-of-plane directions at 
sufficiently small temperatures, protected from switch- 
ing to the in-plane perpendicular to the chain direction 
by a value of ss 1.2 meV/Fe. 



The two different magnetic types of chains are very 
close in total energy, which facilitates their experimen- 
tal observation and provides a considerable challenge for 
experimentalists to verify their magnetic ground state. 
With our calculations we provide theoretical evidence 
for the feasibility to use spin-polarized STM to resolve 
the atomic arrangement and magnetic order. Consider- 
ing the rather small calculated total energy differences 
between the different magnetic collinear solutions of the 
order of 20—30 meV/Fe in both types of chains, we also 
cannot exclude the occurrence of noncollinear magnetic 
states either due to exchange interactions or due to the 
Dzyaloshinskii-Moriya interaction driven by spin-orbit 
couplingP^I Such non-collinear calculations were beyond 
the scope of the present work due to the large supercell 
required for realistic modeling of chains on the (5 x 1) 
reconstructed Ir(001) surface. However, future investi- 
gations of this system need to address this open ques- 
tion. Qualitatively and quantitatively different energy 
landscapes of the magnetization direction in real space 
could also result in a different response of the magneti- 
zation with respect to external magnetic field or temper- 
ature, providing an additional channel for tackling the 
magnetism in these two types of chains experimentally. 
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FIG. 7: (color online) Partial charge density plots for the FM ground state of the CI structure (with 7 layers of Ir substrate) 
with an out-of-plane easy magnetization axis in an energy regime of (Ef,Ef + 0.5) eV. (a) shows a cross-section plot of the 
majority states and (b) the minority states up to a distance of 5 A from the Fe chains (c.f. Fig.[l]). (a) and (b) cut through the 
middle green Fe atoms in (f). STM images at a distance of 5 A above the Fe chains are shown in (c) for majority electrons, (d) 
for minority electrons, and in (e) for the sum of both contributions. The part of the surface displayed in the STM images is 
given in (f). Note that the charge density plots are very similar for both spin directions, with the width of the charge density 
around the Fe atoms somewhat smaller for minority electrons. 
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